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1 Introduction 



Until now no experimental or theoretical proof for the absence of the extra 
generations of leptons and quarks has been given. The direct search for extra 
generations and tests of their indirect effects will be one of the major goals 
of the existing and next generation accelerators. For example, CDF [1] gives 
a lower bound of 139 GeV for the fourth generation quark masses assuming 
they are stable enough when leaving the detector. If the interaction with the 
detector is taken into account this lower bound falls by 30-40 GeV. Thus, 
we can barely assume that the lower bound for the masses of the fourth 
generation quarks is somewhere around M z . 

The fourth generation quarks and leptons have been considered in current 
literature due to various theoretical motivations. In [2] a four generation SM 
is proposed to solve a s , R b and R c crises of SM at Z pole. As another solution 
to these problems, in [3] a four generation supersymmetric SM with mt ~ My/ 
is proposed. In addition to these, four generation SM has been used in cal- 
culating the neutron Electric Dipole Moment (EDM) (in the framework of 
the Kobayashi-Maskawa Model) to break the persisting 10~ 32 e — cm theoret- 
ical boundary which is seven orders of magnitude less than the experimental 
upper bound. 

Note that theoretically the smallness of the electroweak radiative correc- 
tions to the LEP1 observables, following Veltman's arguments, enables us to 
conclude that the masses of the fourth generation quarks and leptons must 
be almost degenerate. This is the case because large mass splitting in the 
fourth generation would induce unacceptably large loop corrections to LEP1 
observables. 

In this work we shall base our analysis on a four generation SM. We have 
two basic aims: 1) Determination of those rephasing invariant combinations 
of CKM matrix relevant to the calculation of the neutron EDM. 2) Founda- 
tion of some likelihood bounds for the fourth generation quark masses. 

Denoting the 'top' and 'bottom' quarks of the fourth generation by t' and 
b' respectively, one can find three independent rephasing invariant combina- 
tions of the elements of the CKM matrix: Im[V t * d V tb V* b V cd ], Im[V t * d V tb V t 1 b V t , d ] 
and Im[V t * s VtbV t * b V t ' s ]. Other combinations can be expressed in terms of these 
and the moduli of the elements of the CKM matrix. As it was shown already 
in [4], strange quark dominates in the neutron EDM, and thus a calculation 
of d s yields at least an order of magnitude prediction for the neutron EDM. 
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In this sense our analysis is focused on the calculation of Im[V t *VtbV t * b Vt> s ]. 

For purposes mentioned above, in our analysis we shall make use of the 
experimental and theoretical results on the branching ratios of the FCNC 
decays b — > 57 and b — > sg. 

We will first outline the derivations of the basic expressions for the quan- 
tities of interest. 

Last section is devoted to numerical analysis and discussions. 



2 Derivations 

b — > S7 and b — > sg, being FCNC decays, start occuring at the loop level 
of perturbation theory. These decays involve consecutive b — > u, c, t, t' and 
u, c, t,t' — > s transitions. Thus decay amplitude has the form 

A= \fR{m)/M$,) (1) 

f=u,c,t,t' 

where A/ = Vf s Vf b and the function R represents the combination of the 
results of the loop integrations and the relevant operators responsible for 
b — > 57 and b — > sg decays. The unitarity of the 4x4 CKM matrix gives 

X u + A c + X t + Xtf = (2) 

It is known that A M is much smaller than others, so from (2) it follows that 

A c » -(A t + A t (3) 

Neglecting rri^ and m 2 c compared to M^, combining R's of different flavors 
with the help of (3), and taking the relevant operator structure for the decay 
under concern, one obtains the following expressions for the b — > 37 and 
b — > sg decay amplitudes: 

A b ^ s , = ^=iX t C 7 + X t ,C' 7 ]0 7 (4) 

Ab ^ sg = ^iX t C 8 + X t ,C' s ]0 8 . (5) 

where O7 and 0$ are magnetic and chromo- magnetic penguin operators 

7 = -^m b sa.F P R b (6) 

8 = -^m b scr.G P R b (7) 
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Here Pr = (1 + 7s)/2 and and are photon and gluon field strength 
tensors, respectively. The Wilson coefficients CV, C' 7 , C$ and C 8 in (4) and 
(5) are evaluated at the mass scale /i ~ m b using renorm group equations 
with five active flavours [5,6]. 

CM = -W 6/23 A(x t ) - |(V 4/23 - V 16/23 )B(x t ) + £ (8) 
2 d i=i 



C 8 (/i) = -i( J B(^)- 1.725)r ? 14 / 23 - 0.9135r ? a4086 + 0.0873r ? - - 4230 

-0.057177- - 8994 + 0.020V- 1456 (9) 

where x t = m'/M 2 ,, r, = ^g, and 

hi = (2.2996, -1.088, -3/7, -1/14, -0.6494, -0.038, -0.0186, -O.OO>0) 
at = (14/23, 16/23, 6/23, -12/23, 0.4086, -0.423, -0.8994, 0.1456) (11) 

A(x) and B(x) in (8) and (9) are given by 

x(x -l)(8x 2 + 5x -7) + Qx 2 (2 -3x)ln(x) 
A[X} ~ 12(*-1)4) { ' 

x(x — l)(x 2 — 5x — 2) + 6x 2 ln(x) 
B[X) = 4(x - l) 4 ) ( } 

C7 and C§ could, respectively, be obtained from (8) and (9) by replacing x t 
by xt = m 2 ,/M^. 

In order to minimize the b -quark mass dependence, which leads to 
uncertainities, it is preferable to define braching ratios BR{b — * sj) and 
BR{b — > sg) by BR exp (b — > cez/)/r t?ieor (5 — > cez/) times the decay rates cal- 
culated from the amplitudes in (4) and (5) where Y theor {b — > ceu) is given by 
[6,7] 



r theor( b ^ ceu) = i92 ^ 3 ' ^' ^ mlg(m c /m b )K(m c /m b ) (14) 

with the phase space factor 

g(x) = 1 - 8x 2 + 8x 6 - x 8 - 24x 4 /n(a;). (15) 
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The factor n(x) describes the leading order (LO) QCD corrections and given 
by 

In 

K (x) = 1 - ^((vr 2 - 31/4)(1 - xf + 3/2). (16) 
Finally, upto two -loop accuracy where a s (fj) is given by 

a s{v) = -jrjiX ^27—) (17) 

2 

where L = ln(-df — ), 3 = 23/3 and 3\ = 116/3, for five active flavours. 

QCD 

Then BR(b — > 57) and £>/?(& — > s#) can be written as 

BR(b^ S1 ) = P~f(Cj — C 7 ) 2 [(R + i? 7 ) 2 + J 2 ] (18) 
BR(b^sg) = Pg (C' 8 -C 8 ) 2 [(R + R 8 ) 2 + I 2 ] (19) 

where the various symbols introduced for convenience have the meaning 

p = 6a e I V cs | 2 BR^jp -> cH 
7 ■Kg(m c /m b )K(m c /m b ) 

8a s \V CS \ 2 BR exp (b -> cei/) , , 

ng(m c /m b )K(m c /m b ) 

^ = c^c 7 < 22 > 

H 8 = ™%r (23) 
i2 = Re[-^—] (24) 

At + A t > 

I = Im[-^-] (25) 

Af + A t > 

It is clear that BR(b — > 57) and BR{b — > sg) depend on the combination 
; . Introducing a polar representatinon for A t and \ t > as 

A t = se ie (26) 
A*, = s'e id ' (27) 
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a straightforward calculation yields 



t = / VWTT2 (28) 

S v /(l_ jR )2 + / 2 

Sin(3 = 1 (29) 

,J(R? + P)((1-R)> + P) 

where (3 = 9 — 9'. 

A representation of A t and \ t > as in (26) and (27) is independent of par- 
ticular parametrisation of the CKM matrix [8]. BR(b — > S7) in (18) contains 
theoretical uncertainities coming from V c b, cx s (fj), BR exp (b — > cez/) and quark 
masses. To reduce these uncertainities one can construct the ratio 

_ BRjb -+ sj) 
r ~ B R(b -> s<?) • m 

First solving (18) and (30) simultaneously for R and I and then using 

(28) and (29) one can express s' / s and Sinf3 in terms of r and the parameters 
defined in equatinons (20-23). 

The ratio r in (30) is useful in reducing the theoretical uncertainities. Fur- 
thermore, since (18) and (19) are similar equations in their dependence on 
the parameters with uncertainities (except for a e — > a s interchage) their the- 
oretical calculation in any model cannot produce too different error bounds. 
It is known already that theoretical and experimental results for BR(b — > S7) 
are close to each other assuming specific values for CKM elements [7] . There- 
fore, one can safely regard r to be free of theoretical and experimental un- 
certainities and it can be indentified with the ratio of the central value of 
BR exp (b — > 57) to that of BR(b — > sg). The latter is unknown, except for the 
SM prediction of BR SM (b -> sg) ~ (2.3 ± 0.6) x 10~ 3 [9]. Thus throughout 
the analysis r will be regarded a free parameter. 

^From these observations one can also predict the uncertainities in s' 
and p. ^From (28) one observes that the uncertainity in s' is determined 
mainly by s, namely the uncertainities coming from BR exp (b — > 57) and 
other parameters in (18) contribute negligeably. However, as is seen from 

(29) the uncertainity in (3 is of the same order as that in I. 
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3 Numerical Analysis and Discussions 



Table 1 shows the values of the parameters used in the numerical analysis. 
As is seen from this table BR exp (b — > 57) is uncertain by ~ 30% which 
is much larger than the errors in other quantities. However, as is shown 
in [7] the theoretical result for the branching ratio of the hadronic decay 
BR(B — > X s 7) is uncertain by a similar amount. Thus, as long as BR(B — > 
X s ^) is represented by BR(b — > 57) with the LO Wilson coefficients in the 
framework of the spectator model, the main source of uncertainity comes 
from the experimental result. This causes, in particular, the quantities / and 
R, defined in (24) and (25), be uncertain by a similar amount as BR exp (b — > 
57) . As it was argued at the end of the last section, r in (30) depends mainly 
on the central values of the branching ratios. After analyzing s' and Sin (3 for 
the central values of the quantities in Table 1, the percentage error in their 
values will be set to that of s and /, respectively. 

/,From the Wofenstein parametrisation of the CKM matrix, | X t |~ A 2 , 
where A 0.22 is the Cabibbo angle. In [10] it is argued that Im[\ t \* t ,} ~ A 5 
and \m[V t * d V tb V* b V cd ], lm[V td V tb V t r b V t/d ] ~ A 7 , which makes Im[\ t \* t ,] more 
important than the other two. Thus, to make comparison explicit, in the 
final graphical results we shall present the r dependence of the function g(r) 
defined by 

g(r) = Im[\ t \;}/\ 5 . (31) 

As it will be seen soon, for a fixed value of m t >, Sin/3 is very sensitive to 
the variations in r; Sin/3 is meaningful (real and in [-1, 1] interval) only for a 
limited range of r values. 

If the addition of a fourth generation is an extension of SM, the allowed 
range of r must include the SM prediction of ~ 0.1. This, in particular, 
follows from the sequental character of the families in case of which one does 
not expect the occurence of new operators in the associated decay amplitudes 
as summarized by (4) and (5). In this sense, one can discard some range of 
m t t values that drives the SM prediction outside the window of the allowed 
values of r. 
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3.1 Light Top Case 



The most important non-oblique correction to Z decays arises in the Zbb 
vertex which depends quadratically on m t . Therefore, naively one expects 
Rb to move toward the SM prediction for a lower value of m t . In this context, 
in [3] a supersymmetric SM with m t ~ Mw and m t ' ~ mf DF is proposed. 
Here supersymmetry is necessary to satisfy the requirements of the CDF 
signal. 

In this section we shall analyse g(r) and Sin/? for m t ~ Mw and m ( < ~ 
m t In particular, we shall set m t = 85GeV and vary % in the CDF 
range. This m t and m t i values are away from representing the CDF data. In 
particular, dominance of the decay mode t' — > bW + , and width of t -quark 
must be explained in some appropriate model, such as supersymmetry. We 
leave aside this problem and analyze the light top quark case assuming its 
results might shed light on the construction of some extension of SM. 

In Table 2 we summarize ,as % wanders in the CDF range, the minimum 
ij'min) and maximum (r max ) values of r between which all quantities are 
meaningful. We observe from this table that as increases from 158 GeV 
to 194 GeV the allowed range of r gets narrower. Thus, as m t i approaches to 
the upper CDF limit, r is forced to take values closer to the SM prediction 
of r ~ 0.1. If m t i is further increased r max approaches SM prediction, in 
particular, r ~ 0.1 is just left outside the allowed interval of r values for 
77V ~ 1.57W. 

Fig.l shows the variation of lm[Sin(3] with r for m t = 85 GeV and m t > = 
176 GeV. Sin/5 becomes imaginary for r < 0.027560 and r > 0.16766 leaving 
a window around r ~ 0.1 as the appropriate region. 

Fig. 2 shows the variation of g(r) in (31) with r. It takes the value 3.5 
at r = 0.027560 and vanishes at r = 0.16766. In the rather narrow interval 
0.027560 < r < 0.030 g(r) falls rapidly from 3.5 to unity. Starting from 
unity at r ~ 0.030, it descents gradually to zero at r = 0.16766. For r ~ 0.1, 
g(r) ~ 0.2 making Im[A t A£,] ~ A 6 which is less than the assumtion made in 
[10]. However, for r <~ 0.030, Im[A t A*,] is well above A 5 , approaching A 4 at 
the lower end. 

The errors in the input parameters effect the value of Im[AiAj,] for a 
particular value of r. The allowed range of r , The allowed range of r , 
however, is not sensitive to uncertainities in the parameters in Table 1, as 
it varies with the central value of BR{b — > sg) in a particular model. This 
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analysis is performed for the central values of the input parameters . As 
mentioned before one expects, Sin/5 be uncertain by ~ 30%, so the numerical 
results for Im[A t A£] have an uncertainity of the same order. 

3.2 Heavy Top Case 

In this section we shall analyse g(r) and Sin/3 for m t = 176 GeV and m t i > 
M z . 

In Table 3 we summarize ,as m t / increases from Mz to 300 GeV, the mini- 
mum (r min ) and maximum (r max ) values of r between which all quantities are 
meaningful. We observe from this table that allowed range for r is narrower 
than that in the light top quark case. Infact, as m t i exceeds 200 GeV SM 
prediction for r is just left outside the allowed range of r. As m t ' increases 
further the window of the allowed values of r forgets completely the the SM 
expectation about r. One can therefore fairly say that for heavy top quark 
consistent with CDF results, mass of the top quark of the fourth generation 
must be somewhere in between M z and ~ 200 GeV . 

Fig. 3 shows the variation of lm[Sinf3] with r for m t = 176 GeV and 
m t i = 150 GeV. Sin/3 becomes imaginary for r < 0.02716 and r > 0.11284 
which permits a very narrow interval in which r may take a value. 

Fig. 4 shows the variation of g{r) with r. It takes the value ~ 1.2 at 
r = 0.027116, then falls first rapidly to ~ 0.2 aroud r ~ 0.030, and gradually 
to zero at r ~ 0.11284. As we observe from this figure, g(r) is smaller than 
that in the light top case. Next, one concludes that Im[AtA^,] is less than 
A 5 almost throughout the interval. Especially, around r ~ 0.1, Im[A t Ap] is 
around ~ 0.04, making it ~ A 7 . 

For higher values of m t > the SM prediction of r ~ 0.1 is far outside the 
allowed region. This calculation is performed for the central values of the 
parameters in Table 1. As mentioned before one expects, Sin/3 be uncertain 
by ~ 30%, so the numerical results for Im[AtAj,] has an uncertainity of the 
same order. 

3.3 Comparison of Two Cases and Conclusions 

If one assumes that r does not depart from its SM value, rather clear expec- 
tations follow. In particular, in both cases, it is not possible to draw Im[A t Aj/] 
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to ~ A 5 . As is seen from Fig.2, for r = 0.1, Im[AtA£/] ~ A 6 . In the heavy top 
quark case it is even smaller; as Fig. 4 suggests, for r = 0.1, Im[A t A£/] ~ A 7 ' 4 . 

In both cases, as m t > increases, allowed range for r squeezed to the lower 
end. In particular, for the light top case r = 0.1 point is thrown to "no 
solution" solution for m t > ~ 1.5TeV which is much larger than the m t i ~ 
200Gel / -point of the heavy top quark case. However, choosing large values 
for m t i does not guarantee that g(r) takes values of the order of unity. Except 
for r values around r min , for any m t ', however large, generally g(r) is much 
smaller than unity, forcing Im[A t A*/] << A 5 . Thus, consistency of large m t ' 
and Im[AtAj/] ~ A 5 assumption of [10] is not generally satisfied except for 
the special point of r ~ r m j„. This, in particular, requires a very large 
value for BR(b — > sg) which is unlikely to occur under four -generation SM 
diagrammatics of the b — > sg decay. One thus concludes that, including the 
30% uncertainity in g(r), it is unlikely to have Im[AtA£] ~ A 5 at any m t i in 
both cases. 

In the light top quark case the window of the allowed range of r values is 
wide. Thus, it permits more deviations from the SM expectation. However, 
the SM prediction of r ~ 0.1 is well included in the window. In this sense, 
since r ~ 0.1 is a subset of the allowed r values one has some extension of 
SM. Despite these results, the light top quark case (in this four generation 
SM form) is not acceptable phenomenologically since it is not consistent with 
CDF signal. 

In the heavy top case, the allowed window of r values is narrower than 
that of the light top case. Thus, deviations from the SM prediction is not 
large as in the case of light top case. Width of the window changes rapidly 
with mt'. Moreover, since for m? >~ 200 GeV r ~ 0.1 point is thrown to 
"no solution" region of r values, one can boldly say that m t ' >~ 200 GeV 
is unlikely to be acceptable. It is in this sense that one is able to propose 
some upper bound for m t >. This analysis gives results only on m t /, so fourth 
generation lepton masses and my are not restricted at all. There is no phe- 
nomenological difficulty with this case as long as the "almost- degeneracy" 
condition explained in the introduction is satisfied. If such a fourth sequental 
family of leptons and quarks do indeed exist LEP1.5 or LEP2 will be able to 
detect them. 

One of us (D.A.D.) thanks Scientific and Technical Research Council of 
Turkey for financial support. 
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Parameter 


Range 


^(GeV) 
rrf t DF {GeV) 

1 v cb | 

BR exp (b -> ceu) 
BR ex P(b -> s 7 ) 
A Q CD(GeV) 
M^(GeV) 
m c /m b 


4.8 
176± 18 
0.9743± 0.0007 

(10.4±0.4)% 
(2.32±0.67) 10" 4 
0.195± 0.005 
80.33 
0.30±0.02 



Table 1: Values of the input parameters used in the numerical analysis 



mt>{GeV) 




f'max 


158 


0.027523 


0.18000 


176 


0.027560 


0.16766 


194 


0.027590 


0.15765 



Table 2: Boundaries of the allowed range of r values for m t = 85GeV . 
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m t i{GeV) 




f'raax 


91 


0.0274400 


0.15903 


150 


0.0271600 


0.11284 


200 


0.0269943 


0.09644 


300 


0.0269500 


0.08080 



Table 3: Boundaries of the allowed range of r values for m t = 17QGeV. 
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